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Abstract 

The presence of a nonconformally invariant term in the photon sector of the Lorentz-violating 
extension of the standard model of particle physics, the "Kostelecky term" Ck oc {kF)ai3ij.vF°'^ F'^^ , 
enables a superadiabatic amplification of magnetic vacuum fluctuations during de Sitter inflation. 
For a particular form of the external tensor kp that parametrizes Lorentz violation, the generated 
field possesses a planar symmetry at large cosmological scales and can have today an intensity 
of order of nanogauss for a wide range of values of parameters defining inflation. This peculiar 
magnetic field could account for the presently observed galactic magnetic fields and induces a small 
anisotropization of the Universe at cosmological scales. The resulting Bianchi I model could explain 
the presumedly low quadrupole power in the cosmic microwave background radiation. 

PACS numbers: 98.62.En, ll.30.Cp, 98.80.-k, 98.70.Vc 
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I. INTRODUCTION 

The standard model of particle physics together with 
the Einstein theory of gravity encompasses, basically, all 
known fundamental physics. However, due to the classi- 
cal character of general relativity, these theories are be- 
lieved to be incomplete, and a search for a more fun- 
damental theory which overcomes the stumbling block of 
quantizing gravity is a major goal of present-day theoreti- 
cal physics. A promising candidate is string theory, which 
deals with gravity in a quantum and self-consistent way. 
Although string theory is well away from being an exper- 
imentally tested theory, some of its peculiar, low-energy 
manifestations, if ever detected in astrophysical and/or 
ground-based experiments, could be a strong signal of its 
correctness. Indeed, as pointed out by Kostelecky 
a possible detectable signature of string theory at low 
energies is violation of Lorentz symmetry Q . 

If on the one hand, terrestrial and astrophysical ex- 
periments have not yet either confirmed or ruled out 
the existence of effects of Lorentz violation (LV), on 
the other hand astrophysical observations have definitely 
confirmed the presence of large-scale correlated, micro- 
gauss magnetic fields in any type of galaxies and clusters 
of galaxies. (For reviews on cosmic magnetic fields see 
Ref. Q.) This peculiar property of possessing, roughly 
speaking, the same intensity and correlation scale every- 
where in the present Universe, has been interpreted as a 
strong hint that cosmic magnetic fields are, indeed, relic 
of the early Universe. Obviously, a damning evidence for 
the primordial origin of cosmic magnetic fields would be 
the detection of their effects on the Cosmic Microwave 
Background (CMB) radiation. Until now, however, ex- 
haustive analyses have only given stringent constraints 
on their properties, without finding any evidence of their 
imprints on the CMB [4]. 

The use of the CMB radiation as a probe of the physics 



of the early Universe has proved to be very fruitful in the 
last years. In particular, the high resolution data of tem- 
perature fluctuations of CMB angular power spectrum, 
provided by the Wilkinson Microwave Anisotropy Probe 
(WMAP) [a, HI have, almost definitively, consecrated the 
so-called "A-dominated cold dark matter" (ACDM) as 
the standard cosmological model of the Universe. 

Nevertheless, the 1-, 3-, and 5-year WMAP data dis- 
play at large angular scales some "anomalous" features, 
the most important ones being the low quadrupole mo- 
ment and the presence of a preferred direction in the 
Universe. 

It is extremely important, however, to stress that, due 
to the importance of residual Galactic foreground emis- 
sion, those anomalies in the CMB anisotropy are still 
subject to an intense debate 0|- 

If there is really a problem with the quadrupole mo- 
ment, then its lowness, indicating a suppression of power 
at cosmological scales, may signal a nontrivial topology 
of the large-scale geometry of the Universe [1]. Indeed, 
several possibilities have been proposed in the recent lit- 
erature to understand that suppression (for other 
large-scale anomalies in the angular distribution of CMB 
see, e.g., Ref. 11 1). Recently enough [ll,[l3|, it has been 
shown that a particular case of the simplest anisotropic 
cosmological model, i.e. the Bianchi I model, could ac- 
count for the smallncss of the quadrupole, without affect- 
ing higher multipoles of the angular power spectrum of 
the temperature anisotropy. Such a proposal of an "ellip- 



soidal universe" has been considered also in Ref. [IJ, ll5| . 

Also, the WMAP data display a particular feature 
which has been deeply investigated in the last years: 
a statistically significant alignment and planarity of the 
quadrupole and octupole modes. This seems to indicate 
the existence of a preferred direction in the Universe, 
which has been named "axis of evil" (AE) [l^. Need- 
less to say, there is no space in the isotropic, standard 
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cosmological model for such a type of features. 

In this paper, we investigate the possibiHty that effects 
of Lorentz violation at inflation could be responsible for 
the creation of large-scale magnetic fields possessing pla- 
nar symmetry. As we will see, these fields can have the 
right intensity and correlation length to explain the exis- 
tence of galactic and extragalactic magnetic fields. More- 
over, because of their peculiar symmetry at cosmological 
scales, they could induce a modification of the (isotropic) 
Robertson- Walker metric, in such a way that the result- 
ing cosmological model is well described by the ellipsoidal 
universe model. This, in turn, could naturally account 
for some peculiar and not-yet-explained anomalous fea- 
tures of CMB radiation discussed above. 

The plan of the paper is as follows. In Sec. 11 we 
discuss the generation at inflation of a plane-symmetric 
cosmic magnetic field in the framework of the Lorentz- 
violation extension of the standard model of particle 
physics. Section III deals with the analysis of CMB 
anisotropies including the asymmetric contribution due 
to the presence of the planar field. Finally, we draw our 
conclusions in Sec. IV. Some technical details are pre- 
sented in the appendixes. 



II. LORENTZ- VIOLATING 
ELECTROMAGNETISM AND PLANAR 
COSMIC MAGNETIC FIELDS 

The large correlation scale of cosmic magnetic fields, 
ranging from ~ lOkpc for magnetic fields in galaxies to 
~ IMpc for those in clusters, and the fact that they are 
found to have approximately the same intensity of a few 
microgauss seems to indicate a common and primordial 
origin, probably to ascribe to some unknown mechanism 
acting during an inflationary epoch of the Universe. If 
one takes into account that the collapse of primordial 
large-scale structures enhances the intensity of any preex- 
isting magnetic field of about a factor 10'^ Q , a primeval 
field with comoving intensity of order of nanogauss and 
correlated on mcgaparsec scales could explain the "mag- 
netization of the Universe" . 

During inflation all fields are quantum mechanically 
excited. Because the wavelength A associated to a given 
fluctuation grows faster than the horizon, there will be a 
time, say ii, when this mode crosses outside the horizon 
itself. After that, this fluctuation cannot collapse back 
into the vacuum being not causally self-correlated, and 
then "survives" as a classical real object jl^ . 

The electromagnetic energy density at the time of 
crossing is then fixed by the Gibbons-Hawking tempera- 
ture Tgh 



ter inflation). Taking into account the expression for 
the electromagnetic energy in standard Maxwell elec- 
tromagnetism, one arrives to the result that the spec- 
trum of magnetic fluctuations at the time of horizon 



crossing is given by Bi ~ TJ^ - M^/to|,i [17|, [ij, [l 



rjii 



(1) 



where H is the Hubble parameter (in this paper we con- 
sider, for the sake of simplicity, just the case of de Sit- 



where in the last equality we used the Friedmann equa- 
tion = {8n/3)M'^/mli. Here, is the total energy 
density during inflation (which is constant during de Sit- 
ter inflation) and mp\ ~ lO^^GeV is the Planck mass. 
Because of conformal invariance of Maxwell electromag- 
netism one finds, however, that the present magnitude 
of the infiation-produced field at the scale, say lOkpc, is 
vanishingly small, Bq 10~^^G [l^. (This is true only if 
the background metric is spatially-flat which is the 
case discussed in this paper.) Since the pioneer work of 
Turner and Widrow [l9|, a plethora of mechanisms has 
been proposed for generating cosmic magnetic fields in 
the early Universe, all of which repose on the breaking 
of conformal invariance of standard electrodynamics (see 
references in Ref. ^] and, for recent papers, Ref. [2lL|22]). 

In particular, Kostelecky, Potting and Samuel [23| first 
pointed out that the breaking of conformal invariance is 
a natural consequence of LV. Indeed, they argued that 
the appearance of an effective photon mass, owing to 
spontaneous breaking of Lorentz invariance, could en- 
able the generation of large-scale magnetic fields within 
inflationary scenarios. The idea that Lorentz symmetry 
breaking could result in the generation of cosmic mag- 
netic fields has been pursued since then by others au- 



thors [2J,|2£ 

The aim of this paper is to show that within a par- 
ticular Lorentz-violating model of particle physics, it is 
possible to generate magnetic fields of cosmological type 
possessing a peculiar spatial geometry. Then, in the next 
section, we will analyze their impact on the isotropy of 
the Universe and, in particular, on the cosmic microwave 
background radiation. 

The model we are going to study is the so-called stan- 
dard model extension (SME) [2^, which is an effective 
field theory including all admissible Lorentz-violating 
terms in the Glashow-Weinberg-Salam gauge theory In 
curved spacetimes, the SME action for the photon field, 
here referred to as the Maxwell-Kostelecky (MK) action, 
reads [30|] 

Smk = Jd'^xel-l F^.F^^" - i (fcF)a/3M.i^"^F^1 , (2) 

where F^^, = d^Ai, — d^A^ is the electromagnetic field 
strength tensor and e the determinant of the vierbein. 
The presence of the external tensor {kp)cif3fj.i, breaks 
(particle) Lorentz invariance [30| and parametrizes then 
Lorentz violation. The external tensor {kp)af3^u is fixed 
in a given system of coordinates. Going in different sys- 
tems of coordinates will, generally, induces a change of 
the form of (kF)ai3tiv In the following, we assume that 
the form of {kF)ai3fii' refers to a system of coordinates at 
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rest with respect to the cosmic microwave background, 
the so-called "CMB frame." 

It is worth noting that taking {kp)af3^i, as a fixed ten- 
sor corresponds to have an explicit violation of Lorentz 
symmetry. This could introduce in the theory an in- 
stability associated to nonpositivity of the energy. In- 
deed, working in a flat isotropic universe described by a 
Robertson- Walker metric ds^ = a?{drf' — c?x^), where 
0(77) is the expansion parameter and 77 the conformal 
time, and introducing the electric and magnetic fields as 
Foi = —a?Ei and F^- = eijkO^Bk (Latin indices run from 
1 to 3, while Greek ones from to 3), the electromagnetic 
energy density turns out to be (sij 



£ = i(E2 + B^) + {kpUoja-^E.Ej 
1 4 



(3) 



The positivity of the above quadratic form depends 
on the particular form assumed by the fixed tensor 
(fc_F)a/3/i!/, which is frame-dependent. This apparent 
paradox (frame-dependent positivity of the energy) is 
overcome when considering models in which LV is sponta- 
neously broken. In this case, however, action ^ loses its 
character of generality, since the tensor {kF)af3tiL' is now 
regarded as a vacuum expectation value of some tensor 
field with its own dynamics. For this reason, and follow- 
ing Ref. [3l|, we will take {kp)ai3^i, to be a fixed tensor 
but, at the same time, we will impose positivity of the 
energy in the CMB frame. Needless to say, if in another 
system of coordinates the energy is not positive defined, 
this means that the effective theory with explicit Lorentz 
symmetry breaking becomes meaningless in that frame 
and one needs to consider the full theory with sponta- 
neous Lorentz symmetry breaking in order to get physi- 
cally acceptable results. 



The equations of motion follow from action ^ [31 1 



dnia^E,) - eijkdj{a^Bk) + 
drf[2{kF)io{)ja^'^ Ej — {kp)o.ijkejkia^^Bi] + 
dj[2{kF)ijkoa^^Ek - {kF)ijki^kima^^ B„i] = (4) 

and d,{a'^E,) + di[2{kF)iooja~'^Ej - {kF)ozjkejkia~^Bi] = 
0. The Bianchi identities are ^^(a^B) + V x (a^E) = 
and V • B = 0. 

We are interested in the generation and evolution of 
superhorizon magnetic fields, that is to electromagnetic 
modes whose physical wavelength is much greater than 
the Hubble radius H~^, Aphys ^ H~^, where Aphys — aA 
and A is the comoving wavelength. Since arj ^ H^^, 
introducing the comoving wavenumber k ~ Stt/A, the 
above condition reads |A;?7| <C 1. Observing that the first 
Bianchi identity gives on large scales B ~ krjE, where 
B and E stand for the average magnitude of the mag- 
netic and electric field intensities, and assuming that all 
(non-null) components of {kF)ai3nv have approximately 



the same magnitude, we can neglect in Eq. ([4]), on large 
scales, the terms proportional to the magnetic field. Also, 
the next-to-last term is negligible with respect to the first 
one. Therefore, at large scales, Eq. (jl]) reduces to 

d,-,{a'E,) + d.^[2{kFU^ja-''Ej] = 0. (5) 

Assuming that ||(fcF)iOOill ^ a'' we then have 

{kF)ioo]a^^Ej = a, (6) 

where Ci are constants of integration. 

It is plausible to assume that the tensor {kF)ai3tiv is 
a nonincreasing function of time. Indeed, we can take 
the simple form {kF)iooj oc a"^, with p a non-negative 
real number. In this case, it is clear from Eq. ([6]), the 
first Bianchi identity, and the fact that 77 cx a^^ during 
de Sitter inflation, that the average intensity of the mag- 
netic field grows "superadiabatically," B oc 0^+^. [If the 
elements of (fcF)iooi are either all zero or negligibly small 
in the CMB frame, the inflation-produced field is, in gen- 
eral, small-scale correlated and then not astrophysically 
interesting.] 

In Ref. [SjJ, the case was analyzed where {kF)iooj is 
a constant isotropic tensor, {kF)iooj oc Sij, where Sij is 
the Kronecker delta. In this paper, instead, we study 
the case in which {kF)iooj is "maximally" anisotropic, in 
the sense that only one component, say (A:f)3003i is dif- 
ferent from zero. In this case, only the E^ component 
of the electric field is amplified and this in turn means 
that only the B^ and By components of the magnetic 
field grow superadiabatically during inflation. Accord- 
ingly, the amplified magnetic field will possess a planar 
symmetry, whose plane of symmetry is the xy-plane. 

Moreover, for the sake of simplicity, we will concentrate 
only on a single case, that corresponding to p = 2, so that 
we will assume that 



(fcF)jOOj = a '^S.sSjskF, 



(7) 



where fci? is a constant which gives the magnitude of 
Lorentz violation effect at the actual time, a = qq = 1. 
(It is worth noting that the condition 11(^:^)^00^11 ^ a* 
translates into /cp 3> a®.) As we will see, the choice of 
Eq. ([7|) will correspond to the case of inflation-produced, 
scale-invariant, cosmic magnetic fields. ^ In this case, 



^ Up to today there is no experimental evidence that the tensor 
{kp)apfj^^ is time-dependent or even different from zero (see be- 
low). Accordingly, it is worth stressing that our assumption that 
{kp)cii3^v decreases in time should only be considered as a work- 
ing hypothesis. On the other hand, the choice p = 2 follows 
by the observation that the actual cosmic magnetic fields pos- 
sess, approximately, the same intensity on different scales (that of 
galaxies and that of clusters of galaxies), and by interpreting this 
occurrence as a hint that their spectrum is indeed scale-invariant 
(at least on cosmological scales). 
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the electric and magnetic fields scale as Ez cx and 
Bx oc By cx a^, respectively, while E^, Ey, and B^ scale 
adiabatically. 

Before proceeding further, we estimate the spectrum 
of magnetic vacuum fluctuations generated during de 
Sitter inflation in Maxwell-Kostelecky electromagnetism. 
If, as before, we assume that the only nonzero compo- 
nent of {kF)af3tiiy is just given by (fc_F)3003 = kpa^'^, 
the electromagnetic energy density turns out to be £ = 
i(E^+B^)+fcFa~^-E2. Therefore, at the time of crossmg, 
where \kr]\ ~ 1, and since hp ^ a^, we get £ ~ kpUi^Bf, 
where ai — a{ti) and we used the flrst Bianchi identity. 
Here, Bi stands for the average intensity of the magnetic 
field on the plane xy at the time of crossing. Remem- 
bering that £ ~ H'^, we obtain the spectrum of magnetic 
fluctuations when crossing the horizon, 



Bi 



(8) 



It is worth noting that, in order to have positivity of the 
energy, we are forced to assume kp > 0. 

In order to find the actual value of the infiation- 
produced magnetic field, we have to follow its evolu- 
tion during reheating, radiation, and matter eras. Since 
the resulting analysis is very close to that performed in 
Ref. |i51|, we quote just the final result (full details are 
given in Appendix A): 



Bo_ ( M Y ( Tmi 
nG ^ VlO^GeVy l^MeV 



7(n-l)/2 



10 



-37 



i/2 



(9) 



where Tj^h is the reheat temperature (see Appendix A), 
and n takes the values ±1 according to < Trh or 
> Trh- The temperature T*, defined in Eq. (gSD, is 
the temperature below which electric fields are washed 
out by dissipative effects of the primordial plasma (so 
that magnetic fields evolve adiabatically for T < T*). 
The condition T* > Trh means, indeed, that the mag- 
netic field evolves adiabatically from the end of reheat- 
ing until today and is equivalent to having Trh/McV < 
(10-37A:j.)1/7. 

Upper bounds on ||(A:i?)a^^i,|| come from the analy- 
sis of CMB polarization and polarized light of radio- 
galaxies and gamma-ray bursts: they are, respectively, 

(Although the 



-37 



34|. 



[23, 10-^^ [331, and 10" 
former bound is less stringent than the latter ones, it 
covers the whole portion of coefficient space for Lorentz 
violation. The point-source nature of radiogalaxies and 
gamma-ray bursts, instead, allows us to put constraints 
only on limited portions of coefficient space (3^.) 

Now, it is clear from Eq. ([9]) and Fig. 1 that, for a wide 
range of values of parameters defining inflation (i.e., M 
and Trh) and Lorentz symmetry violation (i.e., kp)^ the 
scaling-invariant, present-day magnetic field can be as 
strong as Bq ~ nG, and then could naturally explain the 
presence of galactic and extragalactic magnetic fields. 
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FIG. 1: The plane-symmetric, scale-invariant, inflation- 
produced magnetic field has an actual intensity of order of 
nanogauss if the values of parameter defining infiation, i.e. 
the energy scale of inflation M and the reheat temperature 
Trh, stay on the curves, kp ~ ||(feF)a/3,ti^|| estimates the mag- 
nitude of the (constant) external tensor which parametrizes 
Lorentz violation [see Eqs. ([2]) and ((Tjl]. 



III. ELLIPSOIDAL UNIVERSE AND THE 
CMB QUADRUPOLE 

The cosmic magnetic field produced according to the 
mechanism discussed in the previous section possesses 
a planar symmetry. The energy-momentum tensor for 
a plane-symmetric magnetic field of the form B = 
{Bx,By,Q) is 
















h{Bl - Bl) 


BxBy 








BxBy 


\{Bl - Bl) 
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-iB^; 



Bj' y - 



Taking an ensemble average, and assuming statistical 
isotropy in the plane of symmetry, gives {BxBy) ~ 
and (Bl) ~ (B^), so that Eq. (dUj) reduces to (Tb)''^ = 
Pb diag(l, 0, 0, — 1), where ps = (B^)/2 is the average 
magnetic energy density. 

We want now to analyze the large-scale evolution of the 
matter-dominated universe filled with a plane-symmetric 
anisotropic component represented by the planar cosmic 
magnetic field. The total energy-momentum tensor is 
then given by 



T^^ = diag(pB -I- p„i, 0, 0, -ps), 



(11) 



where pm is the matter energy density. We then consider 
a cosmological model with planar symmetry, whose most 
general plane-symmetric line element, compatible with 
Eq. d]), is [35] 



ds^ = dt^ - a^{t){dx^ + dy^) - b^{t) dz^, 



(12) 



5 



where a and b are the scale factors. The metric p2)) 
corresponds to considering the xy-plane as a symmetry 
plane. 

Taking into account Eqs. pT|) and Einstein's 
equations read 



2- - =87rG pB 
a b 



Prn), 



a 
a 



a b 
a b 

2 



0, 



(13) 
(14) 
(15) 



where a dot indicates the derivative with respect to the 
cosmic time. 

Since the conductivity of the primordial plasma is very 
high, we are safe in neglecting interaction effects between 
magnetic field and matter In this case, the magnetic 
component of the energy-momentum tensor is conserved, 
Dpi{TBY^ = 0, so that we have 




PB = 0. 



Let us introduce the so-called "eccentricity" 



1 



id 



(16) 



(17) 



The normalization of the scale factors is such that a (to) = 
b{to) = 1 at the present time to- 

In this paper, we restrict our analysis to the case 
of small eccentricities (that is, we consider the met- 
ric anisotropics as perturbations over the isotropic 
Friedmann- Robertson- Walker background) . In this limit, 
from Eqs. (fT3|) -(fT5 l) . we get the following evolution equa- 
tion for the eccentricity: 

d{ee) 



dt 



3H{ee) = SnGpn. 



(18) 



Here H — d/a is the usual Hubble expansion parameter 
for the isotropic universe. In the matter-dominated era, 
it results in a{t) oc so that H — 2/{2>t). Moreover, to 
the zero-order in the eccentricity, from Eq. (fT6|) it follows 
that the magnetic energy density scales in time as cx 



The solution of Eq. UHl) is 



(0) 



1 - 



,3/2 



(19) 



where f2^' — ps(to)/Pcr\ Rev = 'SHq/SttG is the actual 
critical energy density, and Hq = 100 /ikmsec~^Mpc^^ 
is the current Hubble parameter with /i ~ 0.72 the little- 
h constant Q- At the decoupling, t = tdoc, we have 



-"dec 



^^^''^d^c' where edcc = e(idcc) and Zde 



is the red-shift at decoupling Q. Accordingly, we get 



1091 



ede 



4 X 10 



-3 

nG 



If, for instance, we assume for the present cosmological 
magnetic field strength the estimate Bq ~ nG, which 
is compatible with the constraints derived in Ref. [s^], 
with the presence of galactic magnetic fields 0] , and with 
the results obtained in Sec. II, we get an eccentricity at 
decoupling of about Cdcc — 0.4 x 10^^. 

The presence of a nanogauss, plane-symmetric mag- 
netic field at the surface of last scattering affects the 
temperature anisotropics in such a way, as we are going 
to see, to solve the quadrupole problem. 

Expanding the temperature anisotropy in terms of 
spherical harmonics [38| 



AT{9,< 



(21) 



Z=2 m=-/ 



and introducing the power spectrum 




(22) 



the quadrupole anisotropy is defined by the multipole 
£ = 2 



AT2 



(23) 



where (T) ~ 2.73K is the actual (average) temperature 
of the CMB radiation. The quadrupole problem resides 
in the fact that the observed quadrupole anisotropy, ac- 
cording to 3-year WMAP data (see Table I) and 5-year 
WMAP data (see Table II), is in the range 



(AT; 



2^bs 



and 



(AT2)^, 



(210^ 276) ^K^ 



(213 302) ^K^ 



(24) 



(25) 



respectively, while the expected quadrupole anisotropy 
according the ACDM standard model is 



(AT2)j ~ 1252 /iK^, 



if we take into account the 3-year WMAP data, or 



(AT2)i ~ 1207 pK^, 



(26) 



(27) 



according to 5-year WMAP data. 

If we admit that the large-scale spatial geometry of 
our Universe is plane-symmetric with a small eccentric- 
ity, then we have that the observed CMB anisotropy map 



12, ^ 



is a linear superposition of two contributions 
AT = ATa -I- ATi, where ATa represents the tempera- 
ture fluctuations due to the anisotropic spacetime back- 
ground, while ATi is the standard isotropic fluctuation 

(20) 

caused by the inflation-produced gravitational potential 
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TABLE I: The cleaned maps SILC400, WILC3YR, and 
TCM3YR. Note that the values of a2m in this table corre- 
spond to the values of a2m given in Refs. pol . l4ll . |43 | divided 
by (T) ~ 2.73K. Moreover, the values of Re[a2m] and Im[a2m] 
are in units of 10~®. 



Map 


m 


Re[a2m] 


Im[a2m] 


{AT2f/fiK^ Q/10 ^ 







2.75 


0.00 




SILC400 


1 


-0.56 


1.77 


275.8 6.1 




2 


-6.79 


-6.60 









4.21 


0.00 




WILC3YR 


1 


-0.02 


1.78 


248.8 5.8 




2 


-5.28 


-6.89 









1.22 


0.00 




TCM3YR 


1 


0.10 


1.79 


209.5 5.3 




2 


-5.45 


-6.32 





TABLE IL Quadrupole power obtained from the cleaned 
maps "Hinshaw et al. cut sky", WILC5YR, and HILCM5YR. 
Data from Ref. kS\. 



Map 




Q/IO-*' 


Hinshaw et al. cut sky 


213.4 


5.4 


WILC5YR 


242.7 


5.7 


HILCM5YR 


301.7 


6.4 



at the last scattering surface. As a consequence, we may 
write 



(28) 



We want now to analyze the distortion of the CMB ra- 
diation in a universe with planar symmetry described 
by the metric (jl2p in the small eccentricity approxima- 
tion. The null geodesic equation gives that a photon 
emitted at the last scattering surface having energy E^cc 
reaches the observer with an energy equal to Eo{n) = 
{Eo)il + e^3^n|/2), where (£'o) = -Edec/(1 + ^^dec), and 
n = (711,712,^3) are the direction cosines of the null 
geodesic in the symmetric (Robertson- Walker) metric. 

It is worth mentioning that the above result applies 
to the special case where the normal to the plane of 
symmetry is directed along the z-axis. The general case 
where this normal is directed along an arbitrary direc- 
tion in a coordinate system {xcya, zq) in which the 
a^G?/G-plane is the galactic plane, has been analyzed in 
Ref. [13] ■ Closely following [l^l, we perform a rotation 
R = RxiiS) Rz{(p+n/2) of the coordinate system (x, y, z), 
where Rzi^p + 7r/2) and Rx{'&) are rotations of angles 
(/5 -|- 7r/2 and ?9 about the z- and x-axis, respectively. In 
the new coordinate system the z-axis is directed along the 



direction defined by the polar angles (t?, 1^). Therefore, 
the temperature anisotropy in this new reference system 
is 



ATa _ EoiuA) - (Eo) 



1 



2 2 



(T) - (Eo) 
where ua = {Rn)^ is given by 

TT-A (0,(1)) — cos 6 cos — sin 9 sin cos{(j) — ip) . 



(29) 



(30) 



Equations (l29|) and (l30|) . then, give the general expres- 
sion for the temperature anisotropy induced by a planar 
metric whose normal to the plane of symmetry points in 
the direction {d, ip) in the galactic coordinate system. 

From Eqs. ([21]) and ^0^, it follows that only the 
quadrupole terms {£ — 2) are different from zero: 




dec ' 



-e-*'^sin(2^) e, 
4^e-2.^sin2^el 



2 
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(31) 



Since the temperature anisotropy is a real func- 
tion, we have a/ _„j — (— l)™(a; „i)*. Observing 
that a^_^ = (liria^J* [see Eq. mi we get 
a} _^ = (— l)™(a| ,„)*. Moreover, because the standard 
inflation-produced temperature fluctuations are statisti- 
cally isotropic, we will make the reasonable assumption 
that the a\.^ coefficients are equals up to a phase factor. 
Therefore, we can write [131 



021 = -(4,-1)* = \/ ^ 



(32) 



«22 = («2.-2r = y|e^^^Qi 



where < 0i < 2tt, < 02 < 2tt, and < ^3 < 27r are 
unknown phases, and 



or 



Qi ~ 13.0 X 10" 



Qi ~ 12.7 X 10" 



(33) 



(34) 



according to the 3-year WMAP data or 5-year WMAP 
data, respectively. 

We may fix the unknown direction ("iJ, (p) and the ec- 
centricity by solving Eq. which is (for £ = 2) a 
system of 5 equations containing 5 unknown parameters: 
Cdec, V, (1)2, and ^3. Note that it is always possible to 
choose 020 real, so that = 0. 
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TABLE III: Numerical solutions of Eq. (|28|) obtained by using 
the map SILC400. The values of the angles i9, (/!>2, and 03 
are in degrees. 



edoc/10-' 


■& 




02 


03 


Bo/nG 


0.76 


67.9 


5.8 


31.4 


18.7 


1.88 


0.76 


67.8 


174.4 


33.9 


69.4 


1.89 


0.76 


67.8 


354.3 


121.2 


69.8 


1.89 


0.76 


67.7 


185.5 


121.1 


19.5 


1.89 


0.76 


67.9 


50.1 


121.1 


157.9 


1.88 


0.76 


67.9 


130.0 


138.9 


111.9 


1.88 


0.76 


67.7 


310.3 


131.9 


112.1 


1.89 


0.76 


67.7 


229.6 


131.9 


157.6 


1.90 



TABLE IV: As in Table III, but for the map WILC3YR. 





-d 




02 


03 


Bo/nG 


0.74 


66.4 


4.5 


36.7 


30.8 


1.84 


0.74 


66.3 


175.5 


36.8 


72.3 


1.84 


0.74 


66.3 


355.6 


110.4 


72.2 


1.84 


0.74 


66.3 


184.4 


110.4 


31.1 


1.84 


0.74 


66.5 


57.3 


140.0 


160.5 


1.83 


0.74 


66.5 


122.7 


140.0 


112.0 


1.83 


0.75 


66.2 


303.3 


130.8 


112.5 


1.85 


0.75 


66.2 


236.7 


130.8 


160.0 


1.85 



TABLE V: As in Table III, but for the map TCM3YR. 



fidoc/lO"' 


-d 




02 


03 


Bo/nG 


0.78 


69.2 


0.5 


89.3 


46.2 


1.94 


0.78 


69.2 


179.4 


89.0 


52.5 


1.94 


0.78 


69.2 


359.5 


96.6 


52.3 


1.94 


0.78 


69.2 


180.6 


96.9 


49.5 


1.94 


0.78 


69.3 


49.0 


140.5 


154.5 


1.93 


0.78 


69.3 


131.1 


140.4 


116.4 


1.93 


0.78 


69.1 


311.6 


130.6 


116.6 


1.94 


0.78 


69.1 


228.4 


130.6 


154.3 


1.94 



To solve Eq. ([28)1 for £ = 2, we need the observed val- 
ues of the a2,„'s. We use the cleaned CMB temperature 
fluctuation map of the 3-year WMAP data obtained by 
using an improved internal linear combination method as 
galactic foreground subtraction technique. In particular, 
we adopt the three maps SILC400 [4fl], WILC3YR [iH, 
and TCM3YR 42[ . For completeness, we report in Table 
2m corresponding to these maps. 



I the values of a. 

Numerical solutions of Eq 



80 



60 



40 



20 







• SILC400 
▲ WILC3YR 
■ TCM3YR 
O Axis of Evil 



^« ^ < 



o 



50 100 150 200 250 300 350 

1 

FIG. 2: Numerical solutions of Eq. (|28p obtained by using the 
three maps in Table I. Note that b = 90° —1}, and I = ip, where 
(6, /) are the galactic coordinates. The two open circles at 
(b, I) ~ (20°, 130°) and (6, 1) ~ (10°, 110°) define the direction 
of the axis of evil determined in Ref. [31 for the 5- year WMAP 
temperature sky maps in V-band and W-band, respectively. 



maps, are given in Tables III, IV, and V, respectively. In 
Appendix B, we will show that the system (|28p (for i = 2) 
admits 8 independent solutions. Moreover we observe 
that, for each independent solution (edoc, "i?, V'l '/'2, '/'s) 
shown in the tables, there exists another one given by 
(edoc, TT — "i?, <y5 ± TT, (?!>2, <?!>3), where we take the plus sign if 
(p < TT and the minus sign if > tt. 

Looking at Table I and Eq. ([33[) . we see that the val- 
ues of coefficients [020 1 and [021 1 are much smaller than 
Qi. Moreover, comparing Eq. ([55[) with Table I, and 
Eq. ([M[) with Table II, we deduce, roughly speaking, that 
the value of Qi is about twice the value of Q. Assuming 
Qi > |a2oUa2i|, and Qi ~ 2Q, we will show in Ap- 
pendix B that approximate solutions for Cdcc, f?, and tp 
are 



p2 
■^dec 



15(5-f-3V73) 



24 



Qi 



f ~ arcsm 

2 2 



n/6(3n/73- 5) 



79 



referring to the three 



TT Stt Stt 77r 

From Eqs. ^ and ([23|)-([51|) we get 
edcc ^ 0.8 X 10-2, 
which inserted in Eq. ([20[) gives 
Bo - 2 nG. 



As one can check, the approximate solutions Eqs. (|35|). 
([36[) . and ([37[) are quite close to the numerical values. 

In Fig. 2, we plot the numerical solutions 'd versus ip 
which, defining a symmetry axis, singles out a preferred 



(35) 
(36) 
(37) 

(38) 
(39) 
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direction in the CMB map. We adopt the so-called galac- 
tic coordinates system characterized by the galactic lat- 
itude, b, and galactic longitude, I. In our notation, the 
angle b corresponds to 6 = 90° — i? while / = (p. Equa- 
tion ([55]) gives 



(40) 



I ~ 0°, 45°, 135°, 180°, 225°, 315°, 360°. (41) 

From the above equations and Fig. 2, we see that the 
galactic latitude of the symmetry axis is remarkably in- 
dependent on the adopted CMB temperature fluctuation 
map, while galactic longitude is poorly constrained 



6-20°, 

while Eq. ([37|) can be written as 



It is interesting to observe that, according to Ref. [4J], 
the axis of evil points toward 



(&,/)ae^(20°,130°) 



(42) 



if one takes into account the 5-year WMAP data in the 
V-band, and toward 



(6,0ae=^(10°,110°) 



(43) 



if one uses the W-band 5-year WMAP data, so that it 
seems to be very close to the direction 



(6,0 ~ (20°, 135°) 



(44) 



defined by the symmetry axis in our model. 

It is important to stress, however, that the apparent 
alignment of the above two axes could be just an acciden- 
tal fact. Indeed, the axis of evil is determined by a statis- 
tical correlation (alignment and planarity) between the 
quadrupole and octupole modes, while the presence of a 
planar cosmic magnetic field affects only the quadrupole 
moment in the CMB radiation (this is true to the lowest 
order in the eccentricity parameter considered in this pa- 
per; to higher orders, modifications to the standard value 
of octupole intensity can occur, but they are vanishingly 
small) . 

In any case, it is worth noting that there are already inde- 
pendent indications of a symmetry axis in the large-scale 
geometry of the Universe as, for example, those coming 
from the analysis of spiral galaxies in the Sloan Digital 
Sky Survey [45| . 



the evolution of the Universe, the standard cosmologi- 
cal model. This is the celebrated hot big-bang model, 
"equipped" with inflation, dark energy, and cold dark 
matter. Although these last three ingredients have not 
yet been framed in a definitely theoretical model of par- 
ticle physics, from an operative point of view, this cos- 
mological landscape appears satisfactory. 

However, the ubiquitous existence of large-scale corre- 
lated magnetic fields and some presumed anomaly discov- 
ered in the spectrum of the CMB radiation, together with 
the statement of fact that the evolution of the Universe 
(at least) before the Planck era cannot be understood 
without a self-consistent theory of quantum gravity, jus- 
tify the study of both alternative cosmologies (such as 
anisotropic models of the Universe) and models of parti- 
cle physics beyond the standard model. 

Indeed, in this paper, we analyzed the effects on stan- 
dard cosmology of including in the photon sector of the 
standard model a Lorentz-violating term, the "Kost- 
elecky term" Cr = -l{kF)a(3t,uF°''^ Pt^" . Remarkably, 
we found that its presence is responsible for a superadia- 
batic amplification of magnetic vacuum fluctuations dur- 
ing de Sitter inflation. Moreover, the amplified magnetic 
field possesses a planar symmetry at large cosmological 
scales if the external tensor {kp)ap^^^ which parametrizes 
Lorentz violation, assumes a particular (asymmetric) 
form. This peculiar magnetic field could account for 
the presence galactic magnetic fields and induces a small 
anisotropization of the Universe at cosmological scales 
(which can be described by a Bianchi I model). This, in 
turn, could explain the low quadrupole anomaly in the 
CMB radiation. 

Finally, we would like to stress that anisotropic cosmo- 
logical models, such as that induced by the planar cosmic 
magnetic field discussed in this paper, are expected to in- 
duce a certain amount of polarization in the CMB radia- 
tion [4^. Although it is beyond the aim of this paper to 
discuss this kind of problems, we notice that it has been 
shown 47[ that the 3-year WMAP data on large-scale 
polarization could be in agreement with an anisotropic 
model of a universe of Bianchi I type. An appropriate 
analysis of CMB polarization, resulting from the cosmo- 
logical model presented in this paper, is in progress. 



APPENDIX A 



IV. CONCLUSIONS 

Our knowledge of the Universe has greatly improved 
after the discovery of the cosmic microwave background 
radiation. The observation and the analysis of the tem- 
perature fluctuations in this relic radiation, especially 
in the last decade, has confirmed at a surprising level 
of accuracy the canonical theoretical model describing 



In this appendix, we study the evolution during re- 
heating, radiation and matter eras, of the magnetic field 
produced, according to the mechanism discussed in Sec. 
II, during de Sitter inflation. 

After inflation, the Universe enters into the so-called 
reheating phase, during which the energy of the infla- 
ton is converted into ordinary matter. The reheating 
phase ends at the temperature Trh which is less than 
M and constrained as Trh ^ lO^GeV Moreover, 
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CMB analysis requires M < 10~^mpi [l^l, otherwise it 
would be too much of a gravitational waves relic abun- 
dance, and also one must impose that Trh ^ IGeV, so 
that the predictions of big-bang nucleosynthesis (BBN) 
are not spoiled 3, Sl- 
it is worth noting that the condition fc^? ^ is cer- 
tainly fulfilled during inflation and reheating li kp ^ 



a 
To 
10 



where orh — a(rr{H)- Since orh 



where 



10 ■'GeV is the actual temperature [50|, we have 
< lO-'^^ for IGeV < Trh < lO^GeV. 



-126 



< 



'■RH 



Taking into account that the most stringent bound on 
||(fc_F)a/3pi/|| is 10^^'' (as discussed in Sec. II) we get, 
for a wide range of allowed values of kp, that kp ^ 
during inflation and reheating. This, in turn means, ac- 
cording to Eqs. ([6]), d?]) and the first Bianchi identity, that 
a superadiabatic amplification of magnetic vacuum fiuc- 
tuations during these eras takes place. During reheating, 
in particular, taking into account that rj oc a^/^, we have 

After reheating, the Universe enters the radiation- 
dominated era. In this era, as well as in the sub- 
sequent matter era, the effects of the conducting pri- 
mordial plasma are important when studying the evo- 
lution of a magnetic field. They are taken into ac- 
count by adding to the electromagnetic Lagrangian the 
source term j^A^ [l^. Here, the external current j'^, 
expressed in terms of the electric field, has the form 
= (0, (TcE), where CTc is the conductivity. Plasma ef- 
fects introduce, in the right-hand side of Eq. (|4]), the 
extra term —aada^Ei). In this case, it easy to see that 
modes well outside the horizon (assuming that kp ^ a^) 
evolve as 



E oi exp ^— j dr] a^a^ . 



(45) 



Approximating J drj aJ with r\ aJac and using arj 

H^^, we get 



E (X a exp 



2kpH 



(46) 



In radiation era H ~ T^/mpi [50| and, for temperature 
much greater than the electron mass, the conductivity is 
approximately given by Cc ~ T/a [19[, where a is the 
fine structure constant and T the temperature. Then we 
get 



where 



E (X a exp 



GeV 



10- 



57 \ 1/7 



kp 



(47) 



(48) 



This means that for T > T^, we have E <x (which in 
turn gives B oa since 77 cx a in radiation era) , while for 



T < the electric field is dissipated, so the magnetic 
field evolves adiabatically, B oc a~^. [We have assumed 
that kp ^ from the end of reheating until As is 
easy to verify taking into account Eq. this assump- 
tion is certainly satisfied in our case.] 

Finally, evolving along the lines discussed above the 
inflation-produced magnetic field from the time of hori- 
zon crossing until today, and taking into account Eq. ([8|. 
we easily recover Eq. ([9]). 



APPENDIX B 

In this appendix, we solve the system of Eqs. ((28)l for 
£ = 2. The numerical values of the a2m's are listed in 
Table I, while the a^^'s and a^^^'s are given by Eqs. ([3T|) 
and ([32|l. respectively. 

Since it is always possible to choose 020 real, we 
take — 0. Moreover, the temperature anisotropy is 
real so that we have 02, -m — (— l)™(a2m)*, ^2^-™ — 
and4,_,„ = (-l)'"(aL)*. Therefore, the 
system of Eqs. (pSj) reduces to 



020 



6%/5 



Re[a2i 



Im[a2i 



Re[a22] 



Im[a22] 



[H-3cos(2t?)]e^, 



cosip sin(2i9) e^,. 



Qi, 



(49) 



COS02 Ql, 



(50) 



sinv? sin(2i9)e^g^ ~^ V'^ sin(/)2Qi, 



(51) 



cosi^ sin^i? 



COS(/)3Qi, 



(52) 



sin^aQi, 



(53) 



where Qi is given by Eq. ([33]) . We see that these equa- 
tions form a system of 5 transcendental equations con- 
taining 5 unknown parameters: Cdec, V, 4>2, and ^3. 
Solving Eq. with respect to 'd we get two indepen- 
dent solutions: 



(54) 



where 



^ = 1 arccos ( ' ~ ^ ) . (55) 

2 V 3V37re^^^ ) 

Squaring Eqs. and adding side by side, and 

then solving with respect to Lp, we obtain 8 independent 
solutions: 



(56) 
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where 



1 / a7 ± 13 J + /J^ - 72 
- arccos 



and 



7 ^ 




Re[a22] sin^^Je^oc : 
Im[a22] sin^iJe^cc , 



(57) 

(58) 

(59) 
(60) 



By dividing side by side Eqs. ((5T|) and ([50l). and solving 
with respect to 02, we get 



tan(/)2 



/so Im[a2i 



/tF sin^s sin(2'!9) e^^ 



(61) 



/30Re[a2i] + 0Fsin(p sin(2i?) e^^^ 
The same procedure apphed to Eqs. and ([?T|) results 



tan(/)3 



30Im[a22] + a/tt sin(2iy9) sin z^e, 



dec 



30 Re [022] — v^sin(2(p) sin ?9e 



(62) 



dec 



Finally, by squaring Eqs. (|50| and ([STjl . and adding side 
by side, we get 



^dcc ~ 2ce(j(,j. — 0, 
where we have defined 
30 



(63) 



c{ip,'d) = \j — (Re[a2i] cos(^ — Ini[a2i] sin 1^9) csc(2z9), 



diip,{f) = 10[ Qf --\a2i\']csci2^) 



(64) 
(65) 



We observe that the couple {{), ip) can assume 16 different 
values, according to Eqs. ([54|) and ([56]) . Inserting these 
values in Eqs. (|63p - ([S51) we arrive at 16 different equa- 
tions for edcc- It is straightforward to verify that only 
8 of these are "independent," in the sense that, given a 
solution (cdcc,'?, v) of one of the independent equations, 
then (cdcc, TT— t?, (fiztn) is a solutions of one of the "depen- 
dent" ones (we must take the plus sign ii ip < tt and the 
minus sign if (p > n). The 8 independent equations can 
be solved numerically and their solutions are presented 
in Table III, IV, and V. 

Finally, we derive the approximate solutions (|35p -p6 |) . 
To this end, we may formally solve Eq. ([55]) to get 

'^dec 



= c± Vc^ + d. If Qi > |a2i| then d > c^, and 



we obtain 



"'dec 



lOQi |csc(27?)|. 



Inserting the above expression in Eqs. (|54|) - ([55)) . assum- 
ing that Qi ^ 1 020 1, and after same algebraic manipula- 
tion, we get 



|sin(2?9)| 



n/6(3V73- 5) 



79 



(67) 



Inserting the above equation in Eq. ((66|) gives Eq. ([35|) . 
while its solution is indeed Eq. ([36| . 

Assuming that the main contribution to Q comes from 
the 022, and that Re[a22] — Im[a22], we can write 



Re[a22] ^ Im[a22] ^ ~\l j^Q- 



(68) 

Consequently, inserting Eqs. ((66|) - ((68)) in Eq. ([57)) . we get 
1 



(69) 



where 



30(V73- 5) 



0.15 



5760 
2i , 



(H-5\/73)^ -f 120 



Q 



Q 

'Qi 



0.40 



Qi 



(70) 



Observing that the value of the quadrupole according to 
the ACDM standard model is, approximately, twice the 
value of the observed one, Qi ~ 2Q [compare Eq. p3|) 
with Table I, and Eq. ([Ml) with Table II], we roughly 
get (7 ~ 1/2 and, accordingly, Lp+ ~ and ip- ~ 7r/4. 
Inserting these values in Eq. ([55]) . we recover Eq. ([57]) . 
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